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Abstr act. We consider in this paper the rigorous justification of the Zakharov- 
Kuznetsov equation from the Euler-Poisson system for uniformly magnetized 
plasmas. We first provide a proof of the local well-posedness of the Cauchy 
problem for the aforementioned system in dimensions two and three. Then we 
prove that the long-wave small-amplitude limit is described by the Zakharov- 
Kuznetsov equation. This is done first in the case of cold plasma; we then 
show how to extend this result in presence of the isothermal pressure term 
with uniform estimates when this latter goes to zero. 

1. Introduction 
1.1. General Setting. The Zakharov-Kuznetsov equation 

(1) ut + u dr^u + dj:Au = 0, u = u{x, y, z, t), (x, y, z) e W\ t £ R, d = 2,3 

was introduced as an asymptotic model in [28] (see also [17], [15], [6], and [27] for 
some generalizations) to describe the propagation of nonlinear ionic-sonic waves in 
a magnetized plasma. 

The Zakharov-Kuznetsov is a natural multi-dimensional extension of the Korteweg 
de Vries equation, quite different from the well-known Kadomtsev-Petviashvili (KP) 
equation though. 

Contrary to the Korteweg-de Vries or the Kadomtsev-Petviashvili equations, the 
Zakharov-Kuznetsov equation is not completely integrable but it has a hamiltonian 
structure and possesses two invariants, namely (for uq — u(-,0)) : 

(2) M{t)= I u^{x,t)= I ul{x)^M{0) 
and the hamiltonian 

(3) ^W-J/ [|Vu|^-t] = ^/ [|VuoP-f] = ^(o)- 

The Cauchy problem for the Zakharov-Kuznetsov equation has been proven to 
be globally well posed in the two-dimensional case for data in 7J^(M^) ([7j), and 
locally well-posed in the three-dimensional case for data in H^iJA^), s > |, ([H]) 
and recently in //"(K^), s > 1, [21] ■ We also refer to [22] for solutions on a nontrivial 
background and to [20] and [25] for well-posedness results of the Cauchy problem for 
generalized Zakharov-Kuznetsov equations in R^. Unique continuation properties 
for the Zakharov-Kuznetsov equation were established in [33] , [2] . 
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The Zakharov-Kuznctsov equation was formally derived in |28| as a long wave 
small-amplitude limit of the following Euler-Poisson system in the "cold plasma" 
approximation, 

dtn + V - v + V ■ (nv) = 0, 

(4) { 9tv + (v V)v + V(/) + aeAv = 0, 

- e"^ + 1 + n = 0. 

Here n is the deviation of the ion density from 1, v is the ion velocity, </> the electric 
potential, a is a measure of the uniform magnetic field, applied along the vector 
e = (1, 0, 0)-^ so that if v = (wi, 1^2, ws)"^, e A v = (0, — W3, U2)"^- 
Note that this skew-adjoint term is similar to a Coriolis term in the Euler equations 
for inviscid incompressible fluids. 

The main goal of the present paper is to justify rigorously this formal long- 
wave limit. The one-dimensional case (leading to the Korteweg-de Vries equation) 
has been partially justified in [5], and Guo-Pu gave in a recent preprint a full 
justification of this limit 

Setting p = (1 + n), (jH) writes 

dtp + \7- (pv) = 0, 

(5) { atv + (v- V)v + V(/) + aeAv = 0, 

Ac/) - e"^ + p = 0. 

In this formulation, the Euler-Poisson system possesses a (formally) hamiltonian 
conserved energy 

i/(p,v,0) = / [ip|v|2+p^_i|v</.p_e^ + l]rfa; 

(6) = / [ip|v|2 + i|V0|2+e^0-l) + l]dx, 

as it is easily seen by multiplying the first (resp. second) equation in ([5]) by ^|vp-|-(/) 
(resp. pv), then adding and integrating over M''. 

Of course one has to find a correct functional setting in order to justify the definition 
of H and its conservation (see Remark 1 in Subsection 2.1 below). 
Note that the Euler-Poisson system with a — has another formally conserved 
quantity, namely the impulse 

(7) ^ [ p^dx = 0. 

at j^d 

This conservation law is easily derived (formally) from the equation for pv: 

(8) dt{pv)+V ■{pv(g>v)+pV(f) = 0. 

Linearizing Q around the constant solution n = 1, v = 0, (/) = 0, one finds the 
dispersion relation for a plane wave e'^"*""^-^), k = (ki, k2, fcs): 

(9) ..^(k) + .2(k)(a^ - J^) - . 
or 
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where |kj^p = fc| + k^. 

In the absence of applied magnetic field (a =0), this relation reduces to: 
(11) .\l.)= 1^1' 



l + |k|2- 



Those relations display the weakly dispersive character of the Euler-Poisson system. 

Contrary to the KP case (see 1 for the justification of various asymptotic models 
of surface waves) the rigorous justification of the long-wave limit of the Euler- 
Poisson system has not been carried out so far (see however [H [11] in the one- 
dimensional case). This justification is the main goal of the present paper. 

To start with, we investigate the local well-posedness of the Euler-Poisson system 
dH) which does not raise any particular difficulty but for which, to our knowledge, 
no result seems to be explicitly available in the literature in dimensions 2 and 3. 
The paper [26] concerns a "linearized" version of ([4]), namely the term A0 +1 — e'^ 
is replaced by its linearization at = that is {A — l)(f> (see also [M] where 
a unique continuation property is established for the one- dimensional version of 
this "modified" Euler-Poisson system). Note that this "linearized" version of (|3|) 
is somewhat reminiscent of the pressureless Euler-Poisson system which has been 
intensively studied (see for instance [4], [3] and the references therein). 

In the one-dimensional case (jlj) has the very simple form 

5tn+[(l + n)v], =0, 
(12) { dtv + vv^ (/-x = 0, 

t>xx - e"^ + I + n = 0. 

The existence of supersonic solitary waves for (|T2t has been proven in ^9j. The 
linear stability of those solitary waves was investigated in [13] and their interactions 
was studied in jl2l, in particular in comparison with their approximations in the 
long wave limit by KdV solitary waves. 

Though our analysis is mainly concerned with the higher dimensional case d = 
2,3, our results apply as well to the system (jl2l) and to provide an alternative proof 
to [TT] of the justification of the KdV approximation. 

Let us mention finally, that Q is valid for cold plasmas only; in the general case, 
an isothermal pressure term must be added, 

{9tn + V • V + V • (nv) = 0, 
Vn 
dtv + (v • V)v + V0 + a h ae A V = 0. 
1 + n 
A0 - e"^ + 1 + ?i = 0, 

with a > 0. For this system (with a > 0), global existence for small data has been 
proved [lOj in dimension d = 3 in absence of magnetic field (a = 0) and in the 
irrrotational case. Still for a > 0, in [8], the authors provide for the full equations 
uniform energy estimate in the quasineutral limit (i.e. e'^ = 1 + n) for well prepared 
initial data (note that they also handle the initial boundary value problem). The 
derivation and justification of the KdV approximation is generalized in |11| using a 
different proof as in the case a = 0. We show how to extend our results on the ZK 
approximation to (|13|) and provide uniform estimates with respect to a that allow 
one to handle the convergence of solutions of p3|) to solutions of ([4]) when a 0. 
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1.2. Organization of the paper. The paper is organized as follows. 
In Section [2] we prove that the Cauchy problem for the Euler-Poisson system ^ is 
locally well-posed. The main step is to express (/) as a function of n by solving the 
elliptic equation in (j4]) by the super and sub-solutions method. This step is of course 
trivial when one considers the "linearized" Euler-Poisson system. Then we establish 
the local well-posedness for data (no,vo) in H''-^(R'^) x H^R'^Y, s > | -h 1 such 
that \no\oQ < 1, thus generalizing the result of [26]. We use in a crucial way the 
smoothing property of the map n i— >■ 

In Section [3] we derive rigorously the Zakharov-Kuznetsov equation as a long 
wave limit of the Euler-Poisson system. In order to do this, we need to introduce a 
small parameter e and to establish for a scaled version of the Euler-Poisson system 
existence and bounds on the correct time scale. However the elliptic equation for 
(j) provides a smoothing effect which is not uniform with respect to e and this 
makes the Cauchy problem more delicate (we cannot apply the previous strategy 
which would give an existence time shrinking to zero with e). We are thus led 
to view the Euler-Poisson system as a semilinear perturbation of a symmetrizable 
quasilinear system and we have to find the correct symmetrizer. We obtain in this 
framework well-posedness (with uniform bounds on the correct time scale) for data 
in H'iR'^) X i?|+i(R'')'*, s > f + 1 (we refer to ^ for the definition of H^+^(R'^)). 
We then prove that this solution is well approximated by the solution of the ZK 
equation on the relevant time scales. 

We finally show in Section |4] how to modify the results of Section |3] when the 
isothermal pressure is not neglected, i.e. when one works with ()13|) instead of (j4|. 
The presence of a nonzero coefficient a > in front of the isothermal pressure 
term induces a smoothing effect on the variable n; however, this smoothing effect 
vanishes as a -> 0, and, in order to obtain an existence time uniform with respect 
to e and a, it is necessary to work in a Banach scale indexed by these parameters 
and that is adapted to measure this smoothing effect. 



1.3. Notations. - Partial differentiation are denoted by subscripts, d^, dt, dj = d^j 
etc. 

- We denote by |-|p (1 < p < oo) the standard norm of the Lebesgue spaces 
LP{R'^) (d = 2,3). 

- We use the Fourier multiplier notation: f{D)u is defined as J^{f{D)u){£_) = 
/(^)m(^), where J" and stand for the Fourier transform. 

- The operator A = (1 — A)^/^ is equivalently defined using the Fourier multiplier 
notation to be A = (1 |Z)p)i/2, 

- The standard notation iJ'*(]R''), or simply i?" if the underlying domain is clear 
from the context, is used for the L^-based Sobolev spaces; their norm is written 

- For a given Banach space X we will denote \-\x.t the norm in C{[0,T]; X). 
When X = L^, the corresponding norm will be denoted \-\p,T- 

- We will denote by C various absolute constants. 

- The notation A + {B)s>to refers to ^ if s < io and A + B ii s > tg. 
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2. The Cauchy problem for the Euler-Poisson system 

The aim of this section is to prove the local well-posedness of the Cauchy problem 
associated to the Euler-Poisson system. 

dtn + V -v + V ■ (nv) ^ 0, 

dtv + (v • V)v + V0 + ae A V = 0, 

(14) < A0-p> + l + n = O 

n(-, 0) = no, 
. v(-,0) -vo. 

2.1. Solving the elliptic part. We consider here, for a given n, the elliptic equa- 
tion 

(15) L(0) = -A(/i e"^ - 1 = n. 

Proposition 1. Let n € L°°{R'^) n H^iR'^), d = 1,2, 3, such that infad l+n> 0. 
Then there exists a unique solution </> G i?^(R'^) of US]) such that: 

(i) The following estimate holds 

-K^ = ln(l - \n\oo) <4><K+^ ln(l + \n\^). 

(ii) Defining c.^{n) = |(1 + n)^"'"|oo and Ii{n) = ^ l^li + I'^l/j-i' ''"'^ ^'^^ 

+ 1 |V0P + |A0p + 2e^|V<^p + i(e^ - 1)^] < l/i(n). 

(Hi) If furthermore, uE H''{R'^), s>0, then G H''+'^(R'^) and 

(16) |A^+iV0|2 < Fs{h{n), |nU)(l + Hh^), 
where Fs{., .) is an increasing function of its arguments. 

Remark 1. Writing S,iS,j<j) = —^^\S,\^4' and using the continuity of the Riesz 
transforms we see that we can replace A0 in the left hand side of (ii) in Proposition 
[l]by any derivative d"4>, replacing possibly the right hand side by CIi{n), where 
C is an absolute constant. 



Remark 2. Notice that by (ii), one has I0I2 < 



hjn) 
Coc(n) ' 



Proof. We use the method of sub and super-solutions to construct the solution (p. 
As a super-solution we take 0+ = where is a positive constant satisfying 

K+ > ln(l + |n|oo), 

so that 

£((/)+) > n. 

As a sub-solution, we choose 0_ — —K^ < where is a positive constant 
satisfying 

K- > ln( . / ) = ln(^-), 
mtRd(l + n) Coo(n) 

so that 

L((/«_) < n. 
The next elementary lemma will be useful. 
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Lemma 1. Let he an arbitrary open subset ofW^ and 4> £ L^(r2)ni°°(r2), </> ^ 0. 
Then 

|e^-l|2<^7^ I0I2. 

moo 

If moreover cj) £ H^{n), then e'f' - 1 e H\n). 
Proof. To prove the first point we expand 

n>l n>l 

The last assertion results from the estimate 

\Ve% < el'^l-|V(/)|2. 

□ 

Let F{(j)) — e"^ — 1. Then F'{(p) is bounded by a positive constant K on the 
interval [—K-, K-^-]. Let A > X be so that the function XI — F is strictly increasing 
in [-K-,K+]. 

Let Bp ~ {x e R"^', |x| < p}. We consider the auxiliary problem, for a given cj) e 
H^{Bp) satisfying -K^ <4><K+, 



(17) 



-Alp + Xip^\(j>~ F{cj)) + n 

and we write S{(j)) = ip. Since e L'^{Bp), ip = S{cj)) G H^nH'^{Bp). Moreover, 
one checks easily by the maximum principle that, since —K- — (p- < 4> < = 
<j)+, then also 

-K^ = 0- < 5(0) <K+ = <j)+. 

We now define inductively <f>o — (f)-, (j)k+i = S{(j3k)- By the maximum principle, one 
checks that {(j>k) is increasing, (pk+i > 0fc, A: G N, and that 

0- < 0fc < 0+,A: e N. 

Moreover, S : lF'{Bp) —> L^{Bp) is continuous since z H> — F{z) is Lipschitz. 
The sequence ((pk) is increasing and bounded from above. It converges almost 
everywhere to some cj) which belongs to L'^{Bp) since Bp is bounded. By Lebesgue's 
theorem, the convergence holds also in L^{Bp). On the other hand, using that F 
is Lipschitz on [—K^,K+], one checks that (^fc) is Cauchy in Ill{Bp), proving 
that (j) G Hq{Bp), and is solution of (fT5|) in Bp. Moreover e"^ — 1 € L'^{Bp) and 
(t>£H^{Bp). 

Assuming that (pi and 02 are two Hl\^H^{Bp) solutions we deduce immediately 
that, setting (p = (pi — (p2, 

|V0|^+ / (e*^-e*=)0 = O, 



and we conclude that = by the monotonicity of the exponential. 

To summarize, for any p g N, we have proven the existence of a unique solution 
in Hq n H^{Bp), (which we will denote 4>p from now on), of the elliptic problem 



(18) 



-A^ + e*"- 1 
=0. 
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Moreover, 4>p satisfies the bounds 



(19) 0_ = -K_ <<l>p <(/)+ = K+. 

We derive now a series of (uniform in p) estimates on 0p. We first notice that for 
(j) > —A'-, one has with ao — — , 

(20) (e'^-l)(/.>ao0'. 
Muhiplying by (j>p and integrating over Bp we thus deduce 

(21) 1^ [y\'I^p\' + N^p\']<^HI 

Multiplying ((T8)) by —A<j)p and integrating over Bp we obtain 

(22) [ |A0p|2+ [ e^-lV^pp^ / V0p-Vn. 

Finally we integrate (fT8)) against (e'^p — 1) to get 



(23) / e^^\Vcbp\' + ^ I (e^^-l)2<i 



Adding ([21]), ([221), ([231) we obtain 

[yl^pP + ^|V0pP + |A</.pp + 2e^-|V0pp + \{e^^ - 1)^] 

-Bp 

(24) <l^Hl + \\'^n\l + \\n\l 

Now we extend outside Bp by to get a i/^(M'') function (/)p. Obviously 0p 
satisfies the bound (|21l) and (1231) and up to a subsequence, ^p converges weakly in 
H^iW^), strongly in L^^^(M'^), and almost everywhere to some function (f) S H^{W'') 
which satisfies the bound (IT^ . 

Let us prove that is solution of the elliptic equation L(j) = n (see p5|) ). Let 
X G S(K''). Then supp x'^ Bp for some p. Since (/)p solves L(/)p = n in Bp, one has 



/ Vcj,p-^x+ j {e*' - 1)X = / 



n-X, 

B„ 



and thus 

/ \/(Pp • Vx + / (e"^" - 1)X - / nx- 

jRd jRd JRrf 

One then infers by weak convergence and Lebesgue theorem (using Lemma [T]) that 



/ / (e*-l)x= / nx, 

jRd JK<J jRd 



proving that (j) solves 

(25) - A(/) + e'^ - 1 = n. 

Uniqueness is derived by the same argument used for (jjp. By passing to the limit 
in pT|) and ((23|) one obtains the corresponding estimates for (j>. Since e"^ — 1 — ?i G 
L2(K''), (/) e and ^ for follows. Finally the estimate (ii) in Proposition 

[1] that is (|24p for 0, follows by adding the previous estimates. 
We now prove the higher regularity estimates (iii) assuming that n e i7'*(M'*). 
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From the continuity of the Riesz transforms (see Remark [ij, it is enough to control 
|A'*A(/)|2. We get from ([25]) that 

|A^A0|2 < \n\H^ + - 1\h^ 

the second hne being a consequence of Moser's inequahty. We can therefore deduce 
the result by a simple induction using the first two points of the proposition. □ 

Remark 3. One easily checks that the energy (see ^) makes sense for (rt,v, 0) G 
(R'') xH^{R'^)'^x (R'*)nL°° (R'*) . In fact one has then (recaUing that p = 1+n) , 

P0 - (e^ - 1) = P0 - ((/. + g) = n<t> - .g, <? G L^K^), 
and on the other hand 

p|v|2 = |v|2 + n|v|2 G Li(M'*), 

by Sobolev embedding. 

2.2. Local well-posedness. We establish here the local well-posedness of the 
Cauchy problem for (I14p . 



Theorem 1. Lets > | + no G iJ^^VR"*). vo G H''(M.'^y such thatMl + no > 0. 
There exist T > and a unique solution (n,v) G C([0, T]; iJ''-i(R'') x H^R'^y) of 



([T4| such that 1(1 +n)~^|oo,T < 1 and (j) e C{[0,T];H''+^{R'^)). 
Moreover the energy (j6|) is conserved on [0,T] and so is the impulse (if a — 0). 



Proof. Solving ([T5|) for ip, we set V0 ~ F{n) — VL ^{n) and rewrite (fM)) as 
(26) 



+ V • V + V • (nv) = 0, 
dt-v + (v • V)v + F{n) + ae A V = 0, 
n(-,0) = no, 
v(-,0) = vo. 



We derive first energy estimates. We apply the operator A*^"^ to the first equation 
in (|26|) . A^ to the second one and take the scalar product with A'^^^n and A*v 
respectively. Using the Kato-Ponce commutator estimates 16,^ and integration by 
parts in the first equation, we obtain (note that the skew-adjoint term does not 
play a role here): 



(27) 



i^|A^-in|2 < |A^v|2(l + |A^-in|2) + C|VvU|A^-in|2, 
i||AV|2 < qVvUIA^vl^ + C|A^^^(n)|2|A^v|2. 



As soon as |(1 + n) < oo, using Proposition [T] (for the existence of cj) given 

n) we infer that 

(28) |A'F(n)|2 < Fs-iihin), \n\oo){l + Hh^^i)- 

One has plainly that 

hin) < C\n\l, ( I + 1) < C\n\j,, {- + 1), 
-L — \n\oo Co 
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if 1 + n > f . 



Gathering those inequahties, we obtain the system of differential inequahties for 
\n{-,t)\HB~i and |v(-,t)|/fs, and provided that 1— |n|oo > 



(29) 



1 



I dt 



\n\H,-^ <Ci(|nU.-i,|v|ff.) 
|v|h. <C2(|n|^.-i,|v|ff.) 



where Ci and C2 are smooth functions. Let a{t) — \n{-,t)\Hs-i and /3(t) 
|v(-,i)|_f/s and consider the differential system: 



(30) 



a' = Ci(a,/3)), 
/3' = C2(a,/3). 



Let be the local solution of ([30]) with initial data (Ao,i?o) = {\no\H^-^, \'Vq\h=) 

satisfying 1 + Uq > cq. The solution to (PH)) exists on a time interval which length 
depends only on {Aq,Bq). 

Coming back to (f^ . one deduces that (|«(-,i)|_f/s-i, |v(-,i)|jjs) is bounded from 
above by {A, B) on a time interval / which length depends only on dnol//"-! j |vo|h=) 
(and possibly shortened to ensure that 1 + n > ^ by continuity). 
To complete the proof we need to smooth out This can be done for instance by 
truncating the high frequencies, that is using x(j-D) where x is a cut-off function 
and j = 1, 2, . . . . We obtain an ODE system in x . The energy estimates 

are derived as above and one passes to the limit in a standard way. 

Finally, the conservation of both the energy and the impulse (if a = 0) is obvious 
since the functional setting of Theorem [T] allows to justify their formal derivation. 

□ 



3. The long wave limit of the Euler-Poisson system 

In order to justify the Zakharov-Kuznetsov equation as a long wave limit of 
the Euler-Poisson system with an applied magnetic field, we have to introduce an 
appropriate scaling. 

We set V = [vx^ Vy,Vz)- Laedke and Spatschek [T^ derived formally the Zakharov- 
Kuznetsov equation by looking for approximate solutions of (U) of the form 

0^ = e0(i)(ei/2(^ _ t),e^'^v, e^'^z, e^'H) + e^^^^) + ,3^(3) 
(31) v% = e v(^\e^'\x - t), e^/^y, e'^'z, t^lH) + <?vf> ^vf^ 

v\ = e3/2 4i)(ei/2(^ _ t), ei/2y, e^/^z, <?IH) + e^f + ^5/2^(3)^ 

The asymptotic analysis of the Euler-Poisson system is easier to handle if we work 
with rescaled variables and unknowns adapted to this ansatz. More precisely, if we 
introduce 

X = e"^/^a;, y = e^^^^y, z = e^^^z, i — e^^^t, n — en, (f) — e(j), ir — ev. 
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the Euler-Poisson equation (|4|) becomes (dropping the tilde superscripts) , 

r 9(71 + V • ((1 + e7i)v) = 0, 

(32) I S(V + e(v-V)v + V0 + ae"i/2eAv = O, 
[ -e^A^ + e'-^ - 1 = en, 

and the ZK equation is derived by looking for approximate solutions to this system 
under the form0 

= n(^) {x-t,y,z, et) + eTi^^^ + 6^71^^^ 
= 0(1) {x - t, y, z, et) + e(/.(2) + £20(3) 

(33) vl = vi''' {x - t, y, z, et) + evf^ + e\f^ 

vl = ei/2 „(i) „ + ,^(2) + ,3/2^(3) 

Vl = „(1) (a; _ 2, ri) + + £3/2^(3) 

3.1. The Cauchy problem revisited. It is easily checked that when applied to 
P2p . Theorem [1] provides an existence time which is of order 0(1) with respect 
to e, while the time scale 0(l/e) must be reached to observe the dynamics of the 
Zakharov-Kuznetsov equation that occur along the slow time scale et. 
As explained in the Introduction, we therefore need, in order to justify the Zakharov- 
Kuznetsov equation as a long wave limit of the Euler-Poisson system, to solve the 
Cauchy problem associated to p2l) on a time interval of order 0(l/e). In order to 
do so, we will consider p2p as a perturbation of a hyperbolic quasilinear system 
and give a proof which does not use the smoothing effect of the (j) equation for a 
fixed e. This is the reason why more regularity is required on the initial data in the 
statement of the theorem below. Before stating it, let us introduce the space iJ|+i 
which is the standard Sobolev space H'+'^iW^) endowed with the norm 

(34) Vs>0, V/ei/^+\ = |/|?^.+e|V/|l,.; 

the presence of the small parameter e in front of the second term is here to make 
this norm adapted to measure the smoothing effects of the equation; the fact 
that these smoothing effects are small explains why Theorem [1] which relies on 
them, provides an existence time much too small to observe the dynamics of the 
Zakharov-Kuznetsov equation. 

Theorem 2. Let s > ^ + 1 andna e H''(R'^), vq e H''+^{R'^)'^ such that 1+rio > cq 
for some cq > 0. 

Then there exists T > such that for all e G (0,1), there is a unique solution 
(n%v%0^) e C([0, =];H%'R'^)xH^+^(R'^)'^xH'+\R'^)) of ^ such that I + en" > 
Co/2. 

Moreover the family (ri'^, v*^, V0'^)eg(o.i) is uniformly bounded in H" x _ff|+i x H'^^^ . 

Proof. Step 1. Preliminary results. If we differentiate the third equation of 
([32l) with respect to dj {j = x,y, z,t), we get 

(35) M,{(t))dj(j) = djn with M,(0) = -eA -I- e"^; 

the operator M^{(j)) plays a central role in the energy estimates below, and we state 
here some basic estimates. The first is that {u, M^{(j))u)^^'^ defines a norm which is 



Actually we will not use the third order profiles. 
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uniformly equivalent to more precisely, for all u E H^{M.'^), 

(36) {u,M,{c|))u)<C{\<^\oo)\u\l^ and \u\l^ < C{\cP\^){u, M,{cP)u). 
We also have for all u e H^{R'^) and / e W^-°°iR'^), 

{u, fM,{4>)u) < C(|0U, l/U, Vi|V/|oo) \u\j,^, 

and, for all u G H^{R'^) and / G VT^.^^^j^d)^ 

(u, [/5„M,(0)]u) < C(|</.|h.i,oo, I/Im.1,00, v^|Va,-/|oo) (j = a;,2/,^); 

these two estimates are readily obtained from the definition of Me(<;A) and integra- 
tion by parts. 

Let us finally prove that M^{(j)) is invcrtible and give estimates on its inverse. 

Lemma 2. Let (f) e L°°(R''). Then M^{(j)) : H'^{R'^) -> L'^{R'^) is an isomorphism 
and 

yv e L^iW^), e-5l'^l~|M,((?i)-^t;|2 + \/i|VM,(^)-it;|2 < e5l'^l-|t;|2. 

// moreover to > d/2, s > and (j) e iJ*o+i n H^{R^), then we also have 

e H^R'^), \M,{ct>)-^v\Hs +V^\WM,{cI>)-^v\hs < Cdc^ln^o+in/fOl^^k- 

Proof of the lemma. The invertibility property of M^{(j)) follows classically from 
Lax-Milgram's theorem, and the first estimate of the lemma follows from the coer- 
civity property of M^{4>), namely, 

(37) e-'\^\'-\u\l + e\Wu\l<{M,{<f)u,u). 

In order to prove the higher order estimates, let us write u = M^{<j))~^v. We have 
by definition (— eA + e^'^)u = v, so that applying A'* on both sides, we get 

(38) (-eA + e'"^)k'u = A^w - [A^ e"*']u. 

Using the first estimate, and recalling that u = Mg{(p)^^v, we deduce that 

\M,X^)-\j\hs + Ve\^MM'\\H^ < e'^^^^{\v\Hs + |[A^e'^]A^,(,^)-l^;|2). 

Now, the use Kato-Ponce and Coifman-Meyer commutator estimates, and Moser's 
inequality, shows that 

|[A^e^^]M,(<^)-l^;|2<eC7(|0U) 

X (|?i'|H*o + i|Me(0)-^t;|H,-i + {\(l)\Hs\M^{(l>)~'^v\Hto)s>to + l), 

and we get therefore 

\M,{ct))-'^v\Hs + Vi\VM,{4>)-'^u\Hs < C{\(1)\oo){\v\hs +\(f>\H':o+AMe{(l>)~'^v\Hs-i 

+ {\(j)\Hs\M,{(f>)-\\Hto)s>to+l), 

and the results follows by a continuous induction. □ 

Step 2. estimates for a linearized system. Let T > and (n, v) G 

L°°{[0,T];W^^°°{R'^)) n VFi'°°([0,T];L~(R'^)) be such that 1 + en > cq > uni- 
formly on [0,T], and (p G W^'°° with dt(f) G L°°. Let us consider a classical solution 
(n, V, 0) of the linear system 

{dtn + (1 + en)V • v + ev • Vn = e/, 
dtv + e(v • V)v + V(/> + ae~^/^e A v = eg, 
Me{^)V(l) = Vn + eh. 
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We want to prove here that 

snp {\n\j + < exp(eCor) 

[O.T] 

(40) X (|n|_|^ + Ivi^^Jl^,) + er(|/|i. + |g|^,^ + |h|i. )), 

with Co = Ci^, |(n,v,0)|^yi.oo,|9i(n,v,0)U^,yi|VV • vlioo). 
Taking the L^-scalar product with jq^^, we get 

' -dtfi, n) + (y ■ V, n) + e( — = — • Vti, n) = e( 



1 + en 1 + en 1 + en 

which can be rewritten under the form 

1^/1 \ 1 / dtn , 

7^0t[— n, n) + e-[— ^n, n) + (V • v, n) 

(41) -e^(n, V ■ - )n) ^ e( ^ 

We now take the scalar product of the second equation with Me(0)v to obtain, 
after recalling that M^{(j))W(t) = Vn + eh, 

(M,(0)9tv, v) + e(M,(^)y • Vv, v) + (Vn, v) 
= -e(h,v) +e(Af,(0)g,v), 

or equivalently 

ia*(M,(0)v, v) - ie(a,0e^^v, v) - ie(v, (V • v)Af,(^)v) 

(42) -ie(v, [v • V,M,(^)]v) + (Vn,v) = -e(h, v) + e(M,(^)g, v) . 
Adding ([42]) to (|4T|), we get therefore 



U{-^n,n) + i9,(M,(0)v,v) = + V • (-^)]n,n) 

2 1 + en 2 — 2 (1 + en)^ 1 + en 

+ ie(5t0e^^v, v) + ie(v, (V • v)A/e(^)v) + ie(v, [v • V, A/.(0)]v) 

+e(T^/,^) - e(h,v) +e(A/,(0)g,v). 
1 + en — 

All the terms on the right- hand- side are easily controlled (with the help of Step 1 
for third and fourth terms) to obtain 

x((Y^n,n) + Hl^ + \f\l + + Ih|2)- 

Using ([55)1 and a Gronwall inequality, we readily deduce (HOI) . 

Step 3. H'' estimates for a linearized system. We want to prove here that 

the solution (n, v, (/>) to (|39p satisfies, for all s > + 1, 



sup {\n\l. + |v|^.+0 < exp(eC,r) x |n|^^j2,, + |v|^^J^.+0 

[O.T] ' ^ 

(43) + eT{\f\l^ + |g|^.+, + |h|l,. + \n\l. + |v|^.+0 
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with Cs = C(^, Izilffj, |v|^=+i, |9t(n,v,^)|Loo). 

Applying A'' to the three equations of ([39]). and writing fi — A^n, v — A'^v, and 
(f) = A^cj), we get 

dtfi + (1 + en)W • v + ev • Vn = ef, 
(44) { 5tv + e(v- V)v + V0 + ae"^/^eA V = eg, 

M^{(j))\/4> = Vn + eh. 



with 



/ - A7-[A^n]V•v~[A^v]•Vn, 
g = A^g- [A^v] • Vv, 

h = A"h- i[A",A/,(0)]V0. 
From Step 2, we get that 

sup {\n\l + \^\li) <exp(eCoT) 

[0,T] 

x(l"i..ol2 + \%.o\h) + '^(\f\h + Ig'k. + l^l^i))- 

Now, by standard commutator estimates, we have for ah s > + 1) 

I/I2 + Iglifi < \f\H^ + |g|ff|+i + {\n\H^ + \x\h;+') X (WIh;^' + Hh^) 

and 

|h|2 < \h\H^+Cm„.)\Vct>\H^-i, 

< CmHs){\h\Hs+\n\H^), 

where we used Lemma [2] to get the last inequality. We can now directly deduce 
(|43)) from the Sobolev embedding W^^°°{R'^) C H'iR'^) (s > to + 1). 
Step 4. End of the proof. The exact solution provided by Theorem [1] solves 
([55)1 with (n, V, 0) = (n,v, (/>) and / 0, g = 0, h = 0; we deduce therefore from 
Step 3 that it satisfies the estimate 

sup(|n|l,. + |v|^,+0 <exp(eC,T)x(|n|,^Jl,. + |v|^^J^.+0+^^(l"lff= + l^l?/=+O), 

[0,T] ' ^ ' ^ ' 

with 

Cs = C(— , |v|^=+i, |9t(n,v,0)|Lo?) 

Co 

= C(-^-,|n|ffj,,|v|„a+i), 

Co 

where we used and the equation to control the time derivatives in terms of space 
derivatives, and Proposition [T] to control the L^-norm of (control on V(/) being 
provided by Lemma [5]). This provides us with a T > such that |n||^» + |vP ^+1 
remains uniformly bounded with respect to e on [0,T/e]. □ 
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3.2. The ZK approximation to the Euler-Poisson system. Wc construct here 
an approximate solution to p2p based on tlic ZK equation. Following Laedke and 
Spatschek 19J, but with the rescaled variables and unknowns introduced in (1311) . 
we look for approximate solutions of (|32l) under the form ((33)) . namely, 



= n'-"'"-' (x — t, y, z, et) + en^^-* 

(45) v^^\x-t,y,z,tt) + ev^^'^ 

vl^e^/^v^^\x-t,y,z,et) + evi^l 

Notation 1. We will denote by X the variable x — t and by T the slow time variable 
et. 

Plugging this ansatz into the first equation of ([32l) we obtain: 

6 

(46) dtn^ + W ■ {l + en^)v' =Y^ e^/^N^ , 

3=0 

where 

dx ^ dx^^ 

ATI = + 

dy y dz ' 

dT dx ^ dx^ ^'^dx^ ^ dy y ^ dz ^- 

^3 = |(,(1)41)) + |(„(%(1)) 

iV^ = |;.(^) + A(„(%(2)) + ^(„(l),(2)) 

+ ^("<%i^)) + ^(n(%i^^), 

A.^ = ^(n(%«) + |(n(%«), 

^« = A(.(2),(2))^|(^(.),(.))^|(^(.),(.)). 

Similarly, one has 

o o o ^ ^ 

^ ' dt^ ''dx'' ^dy^ ^dz"" dx^ ^ ^ 



ZK EQUATION 



15 



with Rl =0 and 



for the second component of the velocity equation, we get 



(48) + evl^vl + ev^^vl + e<^«; + - ae'V^,^ = ^e^72ji| 



with 



i!? = -(^4" +««?') 

Rl = i;(2) + A(^;(l)„(2)) +„(1) A^(2) +^(2) 

Rl = t;(2) A^,(2) , „(2) Ai;(2) , „(2) A^;(2) 



while for the third component, the equations are 



6 
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with 



R3-QJ.V. 5^^^. 



i?3 ^ A^,(2) + A„(2) + 1,(1)^^,(1) ^ „(1) Ai;(l) 

«4 ^j,^^. +^ g^^z 



lA„(l) + A/'i;(l)i; f2ll + t,(l)Ai;(2) , „(2) A, 



i?3 _.,(2)_^ (2) , (2)_^ (2) , (2)_^ (2) 

Finally, for the equation on the potential, we obtain 

(50) - e^AcI)' + e"^ - 1 - en = er^ + eV* + eV*^ + ©(e*) 



with 



^2 = ^(1) _ „(1) 



-A0(l)+0(2)^ l(^(l))2_^(2) 
rf^ = -A0(2)+ 0(1)0(2) +i(0(l))3. 

6 

We first derive (following essentially [19]) the equations corresponding to the 
successive cancellation of the leading order remainder terms; we then show that 
they imply that n^^^ must solve the Zakharov-Kuznetsov equation, and then turn to 
show in the spirit of Ben-Youssef and Colin who considered the one-dimensional case 
in [5] (see also jlT]) that it is indeed possible to construct an approximate solution 
(j45p satisfying all the cancellation conditions previously derived. The consequence 
is the consistency property of (|45p stated in Proposition [21 

3.2.1. Cancellation of terms of order zero in e. Canceling the terms N'^ and Rq 
(j = 1,2,3) is equivalent to the following conditions, 

(assuming that <t'''^\ n^'^'^ vanish as |X| +00) 



(51) 






(52) 


"y 






a 


(53) 


z;(i) = 





3.2.2. Cancellation of terms of order e^/^. Using (I5^ - (I5^ . the cancellation of the 
terms and R-{ {j = 1, 2, 3) is equivalent to 

(54) 42)^^5|(i) ^f)^}_Q2 m_ 
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3.2.3. Cancellation of terms of order e. Using the conditions derived above, the 
cancellation of N"^ and {j = 1, 2, 3), 



(55) drn^^^ + 2n^^^dxn^^^ + —dxl^n^^^ = -5x(4'^ - n^^^) 

(56) dx{v^^^ - 0^'') = 9Tn''^ + n(i)9jf 

(57) = 

(58) = 



while the cancellation of r2 is equivalent to (p'^'^^ — n^^^ which has already been 
imposed. 

3.2.4. Cancellation of terms of order e"!"^ . It is possible to cancel the terms of order 
e'^/^ for the equation on the density and on the first component of the velocity; the 
fact that N'^ = B-l = is actually a direct consequence of ([5T |) - ([53l) and ([54| . 
Looking for the cancellation of the other components of the velocity equation, 
namely, setting = = yields respectively 

a 

which are inconsistent with the other equations on n^^^; consequently, we cannot 
expect a better error than 0{e^^'^) on the equations for the transverse components 
of the velocity. 

3.2.5. Cancellation of terms of order . In order to justify the ZK approximation, 
we need to cancel the O(e^) terms in the equation for (j)'^, that is, to impose r"* — 0, 
leading to the equation 

(59) 0(2)-n(2)^An«-^(n«)2. 

3.2.6. Derivation of the Zakharov-Kuznetsov equation. Combining (j55p . (|56l) and 
([59l) . we find that n*^^^ must solve the Zakharov-Kuznetsov equation, 

(60) 2dTn^^'> + 2n'^^^dxn^^'> + + l)Adxn^^'> = 0. 

3.2.7. Construction of the profiles. The ZK equation being locally well posed on 
i/"(M'*), for all s > d/2 + 1,0 we can consider a solution n^^'> e C([0, To]; iJ^(R'^)), 
s > 5 to (|60|) for some Tq > 0. We show here how to construct all the quantities 
involved in (1451) in terms of n^^^ 

• In agreement with ([ST|) . we set (f)^^'^ — vi^"^ = n'^^\ 

• Equations (l52])-(l53l) then give v'y-\vi^^ € CHO^Tq]; H'^-^R'')). 

• We then use §^ to obtain vi^\vP € C{[0,To]; H'-^iW^)). 

• Taking 0(2) ^ J^£>|^„(i) g C([0, Tq]; i?*"2(R'^)) then ensures that ^ 
and (Ull) hold. 

• We then get the density corrector n^^) e C{[0,To]; H''-^(K'^)) by (P)) . 

2 See [T] 1211 124| for the Cauchy theory in larger spaces. 
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• We recover v^^^ £ C{[0,To]; H'-^iR"^)) from (US]) or ([Ml) — the fact that 
n'^^^ solves the ZK equation ensures that we find the same expression with 
(|55l) or ((561) ■ 

The computations above and the expHcit expression of the remaining residual 
terms in (|46|) . (|47| . (|48|) . (j49|) and (|50|) imply the following consistency result. 

Proposition 2. Let Tq > 0, iio e i/"(M'^) (s > 5) and n^^) e C([0, Tq]; i/"(M'^)), 

2aTn(^) + 2n'^^^dxn^^^ +11 + 4) ^dxn'^^'' = 0, n\^'> = uq. 

\ ) 1'=" 

Constructing the other profiles as indicated above, the approximate solution (n*^, v*^, 0*^) 
given by (|45|) solves iS^) up to order in 0*^, m n*^, uj, and up to order c"^"^ in 

r atn' + V • ((1 + en")v') = e^TVS 
(61) i atv' + e(v' • V)v' + V(/)' + ae-i/2eA £3/2^6^ 

i -e^Ac/.' + e'"^' - 1 = en' + e'V^ 

wii/i i?' = (fi^l'^R\,R\,R\) and 

3 

3.3. Justification of the Zakharov-Kuznetsov approximation. We are now 

set to justify the Zakharov-Kuznetsov approximation. 

Theorem 3. Let n° G with s > d/2 + 1, such that 1 + n° > cq on for 

some constant cq > 0. There exists Ti > such that for all e G (0, 1), 

i. The Zakharov-Kuznetsov approximation (ti'jV'j^*^) of Proposition \^ exists on 
the time interval [0,Ti/e]; 

ii. There exists a unique solution {n,v,f) G C{[0,^]; H'iM.'^) x i7|+i(M'')^ x 
H'^'^^{1S.'^)) provided by Theorem\^to the Euler-Poisson equations p2p with initial 
condition (zi", v°. 0°) = (rf 

Moreover, one has the error estimate 

VO < t < Ti/e, \n{t) - n^(t)\l. + |v(t) - v'(t)|^,+, < e^/HC{-,T,, In^U.+a). 

^« Co 

Remark 4. The error of the approximation is 0{c'^'^) for times of order 0(1) but 
of size 0{e^/^) for times of order 0(l/e). Looking at (j45l) . this is a relative error 
of size 0(e^/^) for n and but of size 0(1) for Wj, and w^. Consequently, the 
Zakharov-Kuznetsov approximation provides a good approximation for the density 
and the longitudinal velocity, but not, for large times, for the transverse velocity 
(at least, we did not prove it in Theorem [S]). 

Remark 5. In the one dimensional case (KdV approximation), all the terms of order 
e'^/^ can be cancelled (see §3.2.4|) and the residual in the second equation of (|84|) 
is of size O(e^) instead of 0{e^/^). The error in the theorem then becomes 0{e^t), 
which gives a relative error of size 0(e) for both the density and the velocity over 
large times 0(l/e). In the one dimensional case, it is also possible to construct 
higher order approximations by including the order three and higher terms in the 
ansatz ([33]) . This has been done in pT] . 
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Remark 6. In [11 , the authors justify the KdV approximation (corresponding to 
the one dimensional ZK approximation) by looking at an exact solution as a per- 
turbation of the approximate solution, {uex, Vex) = {napp, Vapp) + ^^{nR, vr), (with 
fc > depending on the order of the approximation) . They study the equations sat- 
isfied by {riR, Vr), which requires subtle estimates. Our approach is much simpler: 
we prove uniform (with respect to e) well posedness of the Euler Poisson equation, 
from which we deduce very easily that any consistent approximation remains close 
to the exact solution (but of course, in a lower norm). 

Proof. Let us take < Ti < min{T, Tq}, where T_/e is the existence time of the ex- 
act solution provided by Theorem[2l and Tp/e the existence time of the approximate 
solution in Proposition H Denote by (n,v,^) e C([0, i7"(M'') x H^+\R'^)'^ x 
Hs+i(^^d-^ the exact solution to ([32]) with the same initial conditions as (n'^, v*^, 0*^) 
furnished by Theorem [51 We also write 

{n,v,(j>) = (n%v%(/>') - {n,v,(j>). 

Taking the difference between (|84p and we get 

{dtn + (1 + en)V • v -|- ev • Vn -t- eVn" • v -|- e(V • v^)n = e/, 
dtv + e(v • V)v + ev • Vv" +W(p + ae^^/^e A v = eg, 
M,{f)\/(j) = Vn - e^e'"^' + eh. 

with / = eiV% g = ei/2i?' and h = eVr' + ~~ ^ - (0 - </)')e"^'') V0". This 

system is of the form p9p with additional linear terms (namely, eVn*^ • v-f e(V • v'^)n 
in the first equation, ev • Vv' in the second one, and —e(j)e^'^ V(/)'^ in the third one) 
that do not affect the derivation of the energy estimate . The only difference 
is that the constant Cs in must also depend on \{Tf ,(j)'^)\^3+i . Since the 
initial conditions for are identically zero, this yields, for all < T < Ti/e, 

sup(|n|2,. -H|v|2 +i) < exp(ea,r) 

[0,T] 

(63) xeT(|/|l,. + |g|^.+, -f |h|l,. -I- |n||,. -I- |v|^.+0' 

with Cs = C(^,|n|Hj,|v|^.+i,Hjy^,|9t(n,v,^)Uoo,|(n%v%0^)|^.+i). Taking if 
necessary a smaller Ti > 0, this implies that for all < t < Ti/e, 

(64) \n{t)\]j. + |v(t)|^.+, < etexp(C.ri) x (|/|^. + |g|^.^, + 
Now, one the one hand, we have 

Cs = C^(-,Nffj,|v|^=+i,H^^j,|9t(n,v,0)U5,,|(n%v%</)^)|^.+i) 

Cq ^' 

= C(-,|nO|^.+5), 
no 

where we used Theorem[2]to control the norms of (n, v, 0) in terms of (n°, v") (with 
v" given in terms of rfi by v" = vj^^_ .^), and the expression of all the components of 

(n*^, v*^, 0*^) in terms of n^^-* to control the norms of the approximate solution. On 
the other hand, we get from the definition of /, g, h and Proposition [2] that 
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We deduce therefore from ([64]) that 

VO<i<ri/e, |7i(t)||,. + |v(t)p <e3/2iC(-,Ti,|nO|^.+5). 

Co 

□ 

4. The Euler-Poisson system with isothermal pressure 

In 11], the authors derived and justified a version of the KdV (and therefore 
d — 1) equation in the case where the isothermal pressure is not neglected. In the 
general of dimension d > 1, the equations ^ are then given by 

{dtn + V -v + V ■ (nv) = 0, 
dtv + (v • V)v + V0 + a-^^ + ae A V = 0, 
A0 - e"^ + 1 + n = 0, 

where a is a positive constant related to the ratio of the ion temperature over the 
ion mass. In the case of cold plasmas considered in the previous sections, one has 
a = 0. In [TT], the cases a — and a > are treated differently, and the limit 
a — > (or, for instance, a — 0{e)) cannot be handled. We show here that the proof 
of Theorem [5] can easily be adapted to the general case a > 0, hereby allowing the 
limit a — > and providing a generalization of the results of [TT] to the case d > 1. 
We first extend Theorem [2] to the general Euler-Poisson system with isothermal 
pressure (|65p . We then indicate how to derive and justify a generalization of the 
Zakharov-Kuznetsov approximation taking into account this new term, in the same 
spirit as the KdV approximation derived in the one-dimensional case in [IL. 

4.1. The Cauchy problem for the Euler-Poisson system with isothermal 
pressure. As in Section [3l we work with rescaled equations. More precisely, we 
perform the same rescaling as for ([32]) : without the "cold plasma " assumption, this 
system must be replaced by 

r 9471 + V • ((1 + e?i)v) = 0, 

(66) I (9tv + e(v ■ V)v + V0 + a + ae'^/^e A v ^ 0, 

I + en 

[ -e^A(l) + e"'' -l^en, 

with a > 0. The presence of the extra term a -^^"^ in the second equation induces 
a smoothing effect that allows the authors of 11 to use the pseudo-differential 
estimates of Grenier [9 . However, these smoothing effects disappear when a ^ 0, 
and the existence time thus obtained is not uniform with respect to a. We provide 
here a generalization of Theorem [5] that gives a uniform existence time with respect 
to e and a (so that solutions to (j66p provided by Theorem |T] can be seen as limits 
when a — ^ of solutions to (l66l) ). This evanescent smoothing effect is taken into 
account by working with n G H^^^(W^) rather than n g H''{M.'^) as in Theorem [2] 
(from the definition ^ of H^^, H^+'^{R'^) coincides with H^R"^) when a = 0). 

Theorem 4. Let s > | + 1, ao > and no e H^+^^{R'^), vq G H''+^{R'^y such 
that 1 — I no loo > Co for some cq > 0. 

Then there exist T > such that for all e G (0,1) and a € (0,Q!o), there is a 
unique solution (n^'",v^^") S C{[0,=]; H^+\R'^) x H^+^{R'^Y) of ® such that 
l + en> co/2 and (j)''" e C{[Q,T]; H'+^{R'^)) . 
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Moreover the family (n*^'", v*^'", V(/)'^'")eg(o,i).Qg(o,ao) is uniformly bounded in H^^^ x 

jjs+l X jjs-l 

Proof. The proof fohows the same steps as the proof of Theorem [2] in addition to 
the operator Me(0) defined in (|35|) . we also need to define another second order 
self-adjoint operator iV^ q,((/), n) as 

(67) 7V,,,(0,n) - -^ + — !— A/,(</.)- ^ 



1 + n 1 + n 1 + n 

provided that infj{£i(l + n) > 0. 

Step 1. Preliminary results. The operator N^^a{<t>Tn) defined in (p7|) inherit 
from the properties of M^{(f>) the following estimates that echoe ((36|) . 



(u,7V,,a((/.,n)w) < C'(|0|oo,|n|oo)|Y- 



en 



1 + en 



^1 <C{\(j)\oo,\Moo){u,Ne^ai4',n-)u)- 



We also have the following commutator estimates, that are similar to those satisfied 
by M^{(f)) (see Step 1 in the proof of Theorem [2]), 

1 I W 1 2 

(u, [dt,N^,a{4','n)]u) < eC( — , \dtn\oo, \n\w^.^ ,\dt(j)\oo, I^U) -— Li 

' ^C(^ ' 1 + en '"'a 

(u, [fdj,N^^a{(l),n)]u) < C( — , |(n,0,/)|vyi,oo, yi^|V9jn|oo,ye|V9j/|oo) 

Co 

I U ,2 



1 + en 



Iff! 



Step 2. estimates for a linearized system. Without the "cold plasma" approx- 
imation, one must replace p9p by 



dtn + (1 + en)V • v + ev • Vn = e/, 

(68) i (9tv + e(v • V)v + V0 + a- h ae^^/^e A v = eg, 

' 1 + en 

M<:(^)V(/) = Vn + eh. 

The presence of the new term in the second equation of (|68p yields some smoothing 
effect on the estimate on n that are absent when a = 0; this smoothing is measured 
with the H^^{M.'^) norm, which coincides with the norm used for (|68l) when 
a = 0. More precisely, we want to prove here that 

sup(|n|2,i +|v||,i) <exp(eCoT) 

[0,T] 

with Co = C ( ^ , I (n , V, 0) I . oo , I a* (n, v, 0) I L~ , v^l v I ^^2, oo , n\^2.^). 

Instead of multiplying the first equation of ((68|) by (1 -I- en)~^ as in the case a — 0, 

we multiply it by N^^aift'Tll) to obtain 

(N,^a{(j>,n)dtn,n) + ([1 + a— ^— M,((/))]V • v, n) + e(iVe,a(</',2i)v • Vn, n) 
— 1 -|- en — — 

= <^{Ne,a{(p,n)f, n). 
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which can be rewritten under the form 
ldt{N,,4^,n)n,n) - U[dt, N,,a{^,n)]n,n) + ([1 + a— ^— Me(^)] V • v, n) 

(70) -e^l'^, [v- V,iV,,a(0,li)]n) = e(iV,,„(0, «)/,«). 

As in the proof of Theorem [2l we take the iP' scalar product of the second equation 
with A/g((/))v; after remarking that 

M,(0)(V(/. + a-^) = [l + aM,(0)-^]Vn + eh 
— ^ 1 + en' — 1 + en 

we obtain with the same computations the foUowing generahzation of ()42p . 
iai(M,(0)v, v) - le(9t^e^^v, v) - ie(v, (V • v)M,(0)v) 

-^e(v, [v • V, Me(^)]v) + ([1 + aMe(0)^-^]Vn, v) 

(71) =-e(h,v) +e(Af,(0)g,v). 

Adding ([71]) to ([70]) . and proceeding exactly as in the proof of Theorem [21 we get 

([Ml). 

Step 3. i/* estimates for a linearized system. We want to prove here that 
for all s > 0, the solution (n, v, (/>) to (|68p satisfies, for all s > + Ij 

sup + |v|^.+i) < exp(eC,T) x (|n|^^J^.+i + |v|^^J^.+i) 



(72) +er(|/|^.+. + |g|^.+. + |h||,. + |n|^.,. + |vp^.,. ) 



with C, ^ C'(^, Inlif^^+i^ Jv|^=+i, , |at(ri, v,0)|hj,). 

Applying A** to the three equations of ([5^ . and writing n = A'^n, v = A^v, and 
(j) = K^(f)^ we get 

dtfi + (1 + en)V • V + ev • Vri = e/, 

(73) <( dtv + e(v • V)v + V0 + a h ae^^/^e A v = eg, 

' 1 + en 

Mt{^)\/(f> = Vn + eh. 

with / and h as in the proof of Theorem [21 while g must be changed into g = 
g + q;[A*, -^_^^^ ]Vn. Using Step 2, one can mimic the proof of the Theorem[2l The 
only new ingredients needed are a control in H^^iW^) of a[A*, -^^^^ ] Vn (the new 

term in g) and a control of / in H^^{M.'^) instead of L^(R''). Classical commutator 
estimates (see for instance [TH]) yield ioi s > d/2 + 1, 

I n I 1 

I/IhL < + {W\m+^ + X (|v|jj.+i + 

so that ([7^ follows exactly as ([^ . 

Step 4. End of the proof. The end of the proof is exactly similar the same as 
for Theorem [2] (the exact solution is no longer furnished by Theorem [T] but by a 
standard iterative scheme). □ 
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Remark 7. The proof given in Theorem [T] for the case a = does not work when 
a > 0. 

4.2. Derivation of a Zakharov-Kuznetsov equation in presence of isother- 
mal pressure. We proceed similarly to the cold plasma case but we replace the 
ansatz P5|) by 

n" = n'-^^ {x - ct, y, z, et) + en'-'^^ 

0^ = 0(1) {x - ct, y, z, et) + e.<lF> 
(74) v'^}\x-ct,y,z,et)^tvf> 

vl^e^'^v^^\x-ct,y,z,et) + ev^^^ 

vl^e^/-'v^^\x-ct,y,z,et) + ev^^\ 
where the velocity c has to be determined. 

Following the strategy of ^13. 21 we plug this ansatz into (|66p , and choose the profiles 
in (|74p in order to cancel the leading order terms. 



4.2.1. Cancellation of terms of order e°. Canceling the leading order 0(1) yields 

-cdxn^^^ + dxv'i^ = 0, 
-cdxv'i'^ + 9x0*^^ + adxn'^^^ = 0, 
(l + a)ayn(i) -aw(i) =0, 
(1 + a)a^n(i) +awW ^ 0. 

4.2.2. Cancellation of terms of order e^/^. We get at this step 



(75) 



(76) 



V, 



(2) ^ 

y 



(l + a 



>3/2 



-dxyH 



(1) 



(1) 



4.2.3. Cancellation of terms of order e. Proceeding as in ii3.2.3[ we get 
(77) dTn^^'> + 2n^^'>dxn'^^'> + \dxAn^^^ 

(2) 



dx {cv 



(2) - 



-9x(4''-cn(2)) 
cdrn'^^^ + n'^^^dxn^^^ 
,(1) 



XXy' 



(l + a)3/%3 



f)i „(i) 



,(1) 



(78) 

(79) 

(80) 9,^ 

(81) ^ 

After replacing (^(^^ by n^^) according to ()8ip . one readily checks that the first two 
equations of ((75|) are consistent if and only if c = \/l + a. 

4.2.4. Cancellation of terms of order e"^"^ . As in § p.2.4p . the cancellation of the 
0(e^/^) terms for the density and the longitudinal velocity equations is automatic, 
but it is not possible for the equations on the transverse velocity. 

4.2.5. Cancellation of terms of order . As in ij3.2.5l we only need to cancel the 
O(e^) terms in the equation for 0', which yields here 



(82) 



0(2)_„(2)^An«-^(n«)2 
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4.2.6. Derivation of the Zakharov-Kuznetsov equation. Combining (|77p. ()78|) and 
((82)) . we find that n*^^' must solve the following Zakharov-Kuznetsov equation 
(which coincides with (1601) when a ~ Q), 

(83) 2cdTn''^^ + 2cn(i)axn(i) + (1 + ^)Mxn^^^ = 0. 

4.2.7. Construction of the profiles. The profiles involved in (1741) are constructed in 
terms of n*^^^ as follows: 

• In agreement with (I5T|) and the first two equations of (175]) . we set ^'^^ = n'^^' 
and v'^^ = cn'^"'^^ with c = + a. 

• The last two equations of (1751 then give Wy^', ui^''. 

• We then use ([7S)) to obtain 

• We take = g2^^^(i) ^ satisfy ([79]) and dSOl). 

• We then get the density corrector n*^^) by ([5^ . 

• We recover from (|77)) or (|75)) - this is equivalent since n'-^^ solves the 
ZK equation (|83)) . 

Finally we get the following consistency result that generalized Proposition [5] 
when isothermal pressure is taken into account. 

Proposition 3. Let Tq > 0, hq e i7"(R'') (s > 5) and n^^^ e C([0, Tq]; i7"(M'')), 
solving ( with c — vT+a ), 

Constructing the other profiles as indicated above, the approximate solution (n*^, v*^, 0*^) 
given by ([7i|) solves \6b]) up to order in (jf , m n*^, v^., and up to order e'^/^ in 

( dtn" + V • ((1 + en')w') = €^N\ 

(84) \ atv' -f e(v^ • V)v^ -t- V# + + ae-i/^e A = e^/^i?^ 

1 + en 

I -e2A^^-)-e^'^■-l = en^ + eV^ 
TOi/i i?^ = (ei/2i?^^,i?|,i?|) and 

3 

i=i 

4.3. Justification of the Zakharov-Kuznetsov approximation. Proceeding 
exactly as for Theorem [3] but replacing Theorem [2] by Theorem SI we get the 
following justification of the Zakharov-Kuznetsov approximation in presence of a 
isothermal pressure. 

Theorem 5. Let n" e with s > d/2 + 1, .such that 1 + en" > cq on R'^ for 

some constant Cq > 0. There exists Ti > such that 

i. The Zakharov-Kuznetsov approximation (n', v', 0*^) of Proposition \^ exists on 
the time interval [0,Ti/e]; 

ii. There exists a unique solution {n,v,cl)) G C{[0,^]; H^+^{R'^) x H^+^{R'^)'^ x 
H''^^{R'^)) provided by Theorem^to the Euler-Poisson equations with isothermal 
pressure (|55)) with initial condition (n^jv",^") = (n^, v', (/)'^)|j^^ . 
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Moreover, one has the error estimate 

Remark 8. The comments made in Remarks 2] and [5] on the precision of the 
Zakharov-Kuznetsov approximation for cold plasmas can be transposed to the more 
general case considered here. 

Remark 9. As already said, the ZK equation (j83p coincides in dimension d = 1 with 
the KdV equation derived in !ll . A consequence of the uniformity of the existence 
time with respect to a in Theorem 2] is that Theorem [5] provides a justification on 
a time scale of order 0(l/e) which is uniform with respect to a whereas it shrinks 
to zero when a — in (llj . 
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